 &MODEL
 NLIST   =           1,
 NSTEP   =        2000,
 EPSI    =  1.0000000E-07,
 WGHT    =   1.000000    ,
 NOBS    =        1000,
 NVAXZ   =           4,
 NVAX    =           3,
 NORDER  =           0,
 NVAZ    =           0,
 NLAMZ   =           0,
 NLAMY   =           1,
 NLAMX   =           3,
 NSCALE  =           1,
 NEXPY   =           0,
 NRES    =           0,
 NCOR    =           0,
 NCOV    =           0,
 VARNAME = BC4      ,
 NVARIANT        =          10
 /
 &YBOXCOX
 VLY     =  0.000000000000000E+000,
 NTLY    =           1,
 NUMLXY  =           0
 /
 &XBOXCOX
 NUMX    =           1,           2,           3, 47*0,
 NUMLX   =           1,           2,           3, 47*0,
 VLX     = 10*0.000000000000000E+000  ,
 NTLX    = 3*1, 7*0
 /
 &YXZDATA
 YXZFILE = TijL2.in                                                                        ,
 YXZFMT  = (5e15.7)                                                                        
 /

VARIABLE       MINIMUM         MEAN          MAXIMUM   NULL OBS.      SCALE
DEPENDENT
Tij         0.5899500D+04  0.2573965D+08  0.3419374D+09     0     0.1000000D-07
INDEPENDENT
Ti          0.1580900D+05  0.1860662D+06  0.3513510D+06     0     0.1000000D-04
Tj          0.6962000D+04  0.1851223D+06  0.4337820D+06     0     0.1000000D-04
Util_BC4    0.3868082D-07  0.7346229D-03  0.1399558D-01     0     0.1000000D+02
CONSTANT    0.1000000D+01  0.1000000D+01  0.1000000D+01     0     0.1000000D+01


MRS of X(l) with respect to X(k) - at the sample means
===================================================================================================================================
        X(k) Ti          Tj          Util_BC4    CONSTANT  
X(l)
===================================================================================================================================

Ti           0.1000D+01  0.1229D+01  0.2776D+09  0.3413D+06

Tj           0.8137D+00  0.1000D+01  0.2259D+09  0.2777D+06

Util_BC4     0.3602D-08  0.4426D-08  0.1000D+01  0.1229D-02

CONSTANT     0.2930D-05  0.3601D-05  0.8135D+03  0.1000D+01


INITIAL ESTIMATES

             LAM    PARAMETER      STD-ERROR        STUDENT-T    GRADIENT       DERIV.Y         ELAST.Y
                                                                                DERIV.E(Y)      ELAST.E(Y)


Ti            X1  0.5325132D+00  0.6598434D-01      8.0702964  0.2799840D-11  0.7366572D+02      0.5325132
                                                                              0.1169046D+00      0.5325132

Tj            X2  0.6510872D+00  0.5906770D-01     11.0227274  0.1252758D-11  0.9052803D+02      0.6510872
                                                                              0.1475469D+00      0.6510872

Util_BC4      X3  0.5837007D+00  0.1616134D-01     36.1171024 -0.3518445D-10  0.2045165D+11      0.5837007
                                                                              0.1703872D+03      0.5837007

CONSTANT      X0  0.9767329D+00  0.1192295D+00      8.1920401  0.7217138D-11  0.2514076D+08      0.9767329
                                                                              0.3300557D+00      0.9767329

VAR(W)            0.1263045D+01  0.5648508D-01     22.3606798  0.0000000D+00

INITIAL VALUES
 
LAM(Y)         0.0000000D+00
LAM(X1)        0.0000000D+00
LAM(X2)        0.0000000D+00
LAM(X3)        0.0000000D+00


LOG-LIKELIHOOD =  0.9344055D+03

NUMBER OF OBSERVATIONS              =           1000

PEARSON-R2(E) (unadjusted)          =         0.2435
PEARSON-R2(E) (adjusted for D.F.)   =         0.2412

PSEUDO-R2(L)  (unadjusted)          =         0.5780
PSEUDO-R2(L)  (adjusted for D.F.)   =         0.5767

PROBABILITY of Y to be at the limit =         0.0000
MEAN VALUE of E(Y)                  =  0.3379181D+00

MEAN of residual W                  =  0.9464949D-14
VARIANCE of residual W              =  0.1263045D+01
95% CONFIDENCE INTERVAL for residual W
Lower Bound = -0.2202751D+01   Upper Bound =  0.2202751D+01




initialization cards for control routine follow - -  id card - no. parameters , no. steps , step size 
minimization tolerance - starting values of parameters - relative uncertainties of parameters (error estimates) -  
and starting direction  (relative parameter increments for 1st step)
identification line
         4      2000   0.1E+01   0.1E-06
  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00
  0.10000E+01  0.10000E+01  0.10000E+01  0.10000E+01
  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00

first entry to fcn , parameters of current step follow - - f,x(i)
  -0.9344054864D+03  0.00000D+00  0.00000D+00  0.00000D+00  0.00000D+00

new absolute minimum , parameters for fcn entry     5  follow - - f,x(i)
  -0.9435799446D+03  0.58114D-01 -0.10201D-02  0.50911D-02 -0.18828D-01

new absolute minimum , parameters for fcn entry     8  follow - - f,x(i)
  -0.9440813130D+03  0.65756D-01 -0.29578D-02  0.15066D-01 -0.32998D-02

new absolute minimum , parameters for fcn entry    13  follow - - f,x(i)
  -0.9604639978D+03  0.61454D-01 -0.29323D+00  0.15312D+01 -0.10163D-01

new absolute minimum , parameters for fcn entry    16  follow - - f,x(i)
  -0.9605955108D+03  0.61610D-01 -0.42753D+00  0.15157D+01 -0.10827D-01

new absolute minimum , parameters for fcn entry    18  follow - - f,x(i)
  -0.9607181967D+03  0.65598D-01 -0.41777D+00  0.15105D+01 -0.16161D-01

new absolute minimum , parameters for fcn entry    20  follow - - f,x(i)
  -0.9607184842D+03  0.65920D-01 -0.41877D+00  0.15120D+01 -0.15927D-01

new absolute minimum , parameters for fcn entry    22  follow - - f,x(i)
  -0.9607184917D+03  0.65916D-01 -0.41878D+00  0.15133D+01 -0.15948D-01

new absolute minimum , parameters for fcn entry    23  follow - - f,x(i)
  -0.9607184917D+03  0.65914D-01 -0.41878D+00  0.15133D+01 -0.15946D-01

new absolute minimum , parameters for fcn entry    24  follow - - f,x(i)
  -0.9607184917D+03  0.65916D-01 -0.41878D+00  0.15133D+01 -0.15948D-01

new absolute minimum , parameters for fcn entry    25  follow - - f,x(i)
  -0.9607184917D+03  0.65915D-01 -0.41878D+00  0.15133D+01 -0.15947D-01


the function (fcn) has not improved by    0.10000D-06 after 3 steps which include a minimum - - -


therefore the run has been terminated
  -0.9607184917D+03  0.65915D-01 -0.41878D+00  0.15133D+01 -0.15947D-01


all numbers of last entry refer to final pt. of run


the lowest value of the function (fcn) was  -0.9607184917D+03 for entry   25 - - x(i) follows
  0.65915E-01 -0.41878E+00  0.15133E+01 -0.15947E-01



run statistics follow - -

 nfcn=   25 nst=   10


MRS of X(l) with respect to X(k) - at the sample means
===================================================================================================================================
        X(k) Ti          Tj          Util_BC4    CONSTANT  
X(l)
===================================================================================================================================

Ti           0.1000D+01  0.2277D+01  0.3688D+09 -0.6420D+08

Tj           0.4392D+00  0.1000D+01  0.1620D+09 -0.2820D+08

Util_BC4     0.2711D-08  0.6173D-08  0.1000D+01 -0.1741D+00

CONSTANT    -0.1558D-07 -0.3546D-07 -0.5745D+01  0.1000D+01


FINAL ESTIMATES

             LAM    PARAMETER      STD-ERROR        STUDENT-T      STUDENT-T    GRADIENT       DERIV.Y         ELAST.Y
                                   (CONDIT.)        (CONDIT.)      (CONDIT.)                   DERIV.E(Y)      ELAST.E(Y)
                                                    (PAR = 0)      (PAR = 1)

Ti            X1  0.1797456D+03  0.6038211D+03      0.2976802                -0.1387997D-06  0.5015131D+02      0.3625327
                                 0.2064242D+02      8.7075850                                0.7671406D+02      0.3948519

Tj            X2  0.2703066D-07  0.1112811D-06      0.2429042                -0.3913793D+01  0.1141822D+03      0.8212107
                                 0.2635639D-08     10.2558253                                0.1171152D+03      0.8901925

Util_BC4      X3  0.1449010D+01  0.4672241D+00      3.1013160                 0.5817389D-06  0.1849765D+11      0.5279324
                                 0.3854853D-01     37.5892352                                0.1893668D+12      0.5500601

CONSTANT      X0 -0.3852034D+03  0.1142939D+04     -0.3370287                -0.5835070D-07 -0.3219766D+10   -125.0897272
                                 0.4889474D+02     -7.8782179                               -0.3282735D+10   -125.8742465

VAR(W)            0.9812504D+01  0.4913118D+01      1.9972050                 0.0000000D+00
                                 0.4194091D+00     23.3960188
 
LAM(Y)    (E)     0.6591481D-01  0.1559440D-01      4.2268254    -59.8987545  0.2332561D-02
LAM(X1)   (E)    -0.4187821D+00  0.2943887D+00     -1.4225481     -4.8194174 -0.4724503D-04
LAM(X2)   (E)     0.1513266D+01  0.3382286D+00      4.4740927      1.5175123 -0.4445941D-04
LAM(X3)   (E)    -0.1594656D-01  0.1561615D-01     -1.0211583    -65.0574246 -0.2871225D-02


LOG-LIKELIHOOD = -0.1745996D+05

NUMBER OF OBSERVATIONS              =           1000

PEARSON-R2(E) (unadjusted)          =         0.2959
PEARSON-R2(E) (adjusted for D.F.)   =         0.2909

PSEUDO-R2(L)  (unadjusted)          =         0.5891
PSEUDO-R2(L)  (adjusted for D.F.)   =         0.5862

PROBABILITY of Y to be at the limit =         0.0000
MEAN VALUE of E(Y)                  =  0.2802748D+08

MEAN of residual W                  = -0.5725654D-09
VARIANCE of residual W              =  0.9812504D+01
95% CONFIDENCE INTERVAL for residual W
Lower Bound = -0.6139684D+01   Upper Bound =  0.6139684D+01


CPU TIME:      0.18 sec
